Decoherence of Rabi oscillations of electronic spin states in a double quantum dot 
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We study the role ol charge fluctuations in the decoherence of Rabi oscillations between spin 
states |tl), |IT) of two electrons in a double dot structure. We consider the effects of fluctuations in 
energy and in the quantum state of the system, both in the classical and quantum limit. The role 
of state fluctuations is shown to be of leading order at sufficiently high temperature, applicable to 
actual experiments. At low temperature the low frequency energy fluctuations are the only dominant 
contribution. 
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The keystone of quantum information processing is the 
coherent dynamics of the quantum logical bits (qubits)^. 
Although such coherent behavior is well established in 
atomic systems, it can be maintained only for very short 
time scales, of the order of few nanoseconds, in charge 
based solid state based systems. To overcome this prob- 
lem one may employ the spin degree of frcdom of the 
electrons residing in a quantum dot as a qubii£. In fact, 
as a consequence of the confined geometry, the coherence 
time of the spin may be extended to be of the order of tens 
of microseconds^, primarily restricted by the coupling of 
the nuclear spin environment via the hyperfine interac- 
tion^. These results motivated the experimental progress 
in controlling electronic spin in GaAs gated quantum dots 
systems^. 

In a recent experiment^ the use of spin states of two 
electrons in a double dot as the holder of quantum infor- 
mation has been investigated. In that configuration the 
system is governed by (i) the hyperfine interaction which 
tends to mix singlet and triplet states and (ii) the ex- 
change interaction which tends to preserve the total spin 
of the electron pair. The interplay between the two ef- 
fects has been studied theoretically^ and analyzed exper- 
imentallyi In particular Ref. [1 reports Rabi oscillations 
between spin states driven (electrostatically) by tuning 
the exchange energy. Such oscillations (faster than the 
typical spin decoherence time) are mainly hindered by 
charge fluctuations^^. 

In this paper we analyze the decoherence effects in the 
Rabi oscillations due to charge fluctuations. We con- 
sider both the effects of exchange energy fluctuations 
and fluctuations of the singlet hybridized state which 
is affected by charge fluctuations as well. In particu- 
lar we calculate the time dependence of the Rabi os- 
cillations in the presence of gate voltage and tunneling 
amplitude fluctuations, both in the classical (high tem- 
perature, Eq. ©) and the quantum case (Eq. jJH]) ). We 
describe the crossover of the decoherence rate between 
low and high temperature regimes, which can be relevant 
in the actual experiments. Classical energy fluctuations 
have been analyzed in Ref. 



The system (Cf. Petta et al£) is schematically pre- 
sented in Fig. [1] It consists of a gate confined semi- 
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Figure 1: (Color online), (a) The schematics of a double 
dot with a nearby QPC; Bnl: Bn r are the respective nu- 
clear magnetic fields, (b) Energy levels of lowest singlet and 
triplet states vs. the detuning parameter e (cf. Ref. 0). (c) 
The time variation of e needed to drive Rabi oscillations be- 
tween |f J.) and HI) states (refer to panel (b)). The system 
is initially in its ground state \S g (eA)}- e is then varied adi- 
abatically (he <C J(0) 2 /Eo) to e = es, keeping the system 
in its lowest S z = state, either ||J.) or \][). During this 
variation the point e = ex is crossed at a time faster than 
H/(gfj,B max{Bjv L , Bnr}) to avoid a transition to |T+) (but 
with he <C J(0) 2 /-Eo); approaching e — cb, instead, the vari- 
ation of e is slowed down (HE i < {T \H N \S g ) 2 < J{0) 2 ) 
to guarantee adiabaticity with respect to the nuclear interac- 
tion. Following the pulse that induces Rabi oscillations in the 
interval [0, r], the adiabatic variation of e at t > r is reversed. 
This allows to identify the state of the system at t = r by 
mapping the states |T!) and |J.f) at e — es to \S g (eA)) and 
| To) at e — ca respectively, where the latter can be measured 
employing the QPC. 

conducting double quantum dot. Tunnel barriers con- 
nect each dot to the adjacent reservoirs allowing dot-lead 
tunneling of electrons. The gate voltages, Vr, Vl, and 
Vr, control the tunnel between the dots, and the dots' 
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charge configuration (tiljUr), respectively. It is possible 
to measure such a charge configuration using a quantum 
point contact (QPC) located near one the dots. The di- 
mensionless detuning parameter, e oc Vl—Vr controls the 
difference nz—UR. In Ref . [H the system was operated be- 
tween (1, 1) and (0,2). In the (0,2) charge configuration 
(e = 6a), the antisymmetric nature of the electron wave 
function enforces the ground state of the system to be a 
singlet. The excitation energy to the lowest triplet state 
is experimentally estimated to be > AQQpeV , larger than 
the charging energy of the single (1, 1) state, E . In the 
(1, 1) configuration (e = er) the singlet and triplet states 
are instead practically degenerate. The low energy states 
of the system are two singlets, \Sl) ,\Sr), correspond- 
ing to (1, 1) and (0, 2) respectively, and the three triplet 
states for charge configuration (1,1), \Tq) , \T+) , |T_), 
respectively with the spin component 0,1,-1 in the z 
direction perpendicular to the dots' plane. An exter- 
nal magnetic field B = Bz is applied to split the states 
\T + ) , |T_) by the Zeeman energy A z = qprB ~ 2.5peV. 
We neglect the role of these states (see below) and we 
write the Hamiltonian of the system in the S z — sub- 
space as 

H = E e(\T )(T \ + \S L )(S L \-\S R )(S R \) 

+E (X s \S L ){S R \+h.c). (1) 

We introduced a tunneling amplitude between the two 
singlet states, A s , which is the only possible tunneling 
matrix element assuming conservation of total spin. It 
can be chosen real and positive. This leads to hybridiza- 
tion of the ground and excited states in the singlet sub- 
space, 

\S g (e)) = - sm9 \S B )+ cos0 \S L ) , 

|S e (e)) = cos9 \S R )+ sm6 \S L ) , (2) 

thence 

Ho = ^(\S e (e)) (S e (e)\-\S e (e)) (S e (e)\)+E e\T ) (T | . 

with ?vjj s = 2E ^/e 2 + X 2 , tan0 = (^e 2 + A 2 + e)/A s . 
The energy levels as a function of e are plotted in Fig. Q] 
Transitions between singlet and triplet states are made 
possible due to the hypcrfinc interaction of the electrons 
in the dot with nuclear spins, which can be written in 
terms of the effective nuclear magnetic field in each of 
the dot, H N = gp B (B NL • Sl + 'Bnr ■ Sr.). The typi- 
cal dynamical scale of the nuclear environment is of the 
order of tens of microseconds, and therefore it acts as if 
it is a frozen external field over the duration of the ex- 
periment. With Bn l, Bn r ~ 5mT -C B, the hypcrfinc 
interaction is effective only at e ~ ex, where it can mix 
\S g ) and and around e ~ cr <C 1 where it mixes 

the lower energy states \S g ) and \Tq): Hn = gpsi^N l~ 
Rnr) ■ z \Tq) (S g \ + h.c, while the energy difference be- 
tween them is J(e) = E {e + ^/e 2 + A 2 ) < (T Q \H N \S g ). 
The ground state of the system at e = er is therefore 



(\S g (e = -1)) ± |T ))/V2 = |U) (for +) or |U). |T!) or 
IT!) ( — ) : tne s P m m the two dots are oppositely oriented. 
Hereafter we consider |||) state. 

In the experiment described in Ref. the detuning 
parameter is varied in time to induce Rabi oscillations 
between |||) and |f|). The time dependence of the pa- 
rameter e used to drive the oscillations is depicted in 
Fig. [TJ The system is prepared in the state |T!) (or 
equivalently |!T)) at e = €r (cf. Fig. [TJ). Subsequently, 
a gate voltage pulse at t = modifies e to a point 
where (To\H^\S g ) <C J(e), thus inducing oscillations be- 
tween |t!) and |!T) over a time interval r. The follow- 
ing manipulation of e (cf. Fig. [T|) allows to relate the 
measured conductance of the QPC with the probabil- 
ity of finding the system in |T!) right after the pulse, 

p ( T ) = I (T!|exp(-rff T) |T!) I 2 - 



III. CLASSICAL NOISE 



The Rabi oscillations are obtained by tuning the en- 
ergy difference J(e) between \S g ) and |Xo), which results 
in the different charge of the triplet and hybridized sin- 
glet. Rabi oscillations will therefore be extremely sensi- 
tive to an environment coupled to charge as opposed to 
the nuclear spin environment. Decohcrence effects will 
originate both from fluctuations of the (exchange) en- 
ergy J(e), analyzed in Ref. [^, and fluctuations of the 
hybridized singlet state \S g (e)). We analyze the Rabi os- 
cillations taking into account fluctuations of Vl , Vr and 
Vt- The respective gates are controlled independently of 
each other, hence it is natural to assume that their fluc- 
tuations are independent. In principle it is possible to 
determine a correlation matrix for the fluctuations of the 
parameters in the Hamiltonian, e and X s , by consider- 
ing a specific potential form for the double dot. Instead 
we assume that Vr affects only the tunneling matrix ele- 
ments, A s , which is reasonable for weak tunneling. Then 
e is affected only by the fluctuations of Vl — Vr. The 
Hamiltonian is 

H = # ({e-e + 6:(t), A,-A,+&(t)}) 

= H Q + E (V c Ut) + VxUt)), (3) 

where are Gaussian distributed with (£i(t)) = 0, 

(ti{t)ZjP)) = 2h 2 T i /Ei5{t-t')6 i , j for i,j G {c,s} and 
V = \T ) (T | + \S L ) (Sl\ - \Sr) (S r \, V x = \S L ) (Sr\ + 
\Sr) (Sl\- The assumption of white noise renders the dy- 
namics Markovian, leading to an exact master equation 
for the density matrix of the double dot, 

d t p = -i/h{H ,p}- £ r j (V J 2 p-2V J pV J + P Vf). (4) 

This differential equation is solved with the initial 
condition p = |T!) (T!| by neglecting terms of order 
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0(J(e)/(hu s ) = sin 6). Under this assumption the den- 
sity matrix at any time t > can be written as 

1 



Pit) = -{\T )(T \+Y(t)\S g )(S g \ + (l-Y(t))\S e )(S e \ 

(5) 



2 Li-"/ \-<ji 1 - vv i-'S/ Wg 

+ (X(t)\S g ){T \+h.c.)} , 



in terms of the two functions Y(t) and X{t) describ- 
ing the evolution of the state populations and co- 
herency respectively. The explicit expressions, X(t) — 
e (iJ(e)/h- j2 )t^ Y (t) = (1 + e" 7lt )/2, allow to deter- 
mine the surviving probability P(t) = p{t) \\\) = 
l/4[l + Y(t) + 2?tte{X(t)}}, 



p ( T ) = g [3 + e^ lT + 4 cos( J^r/^e-^] , (6) 

71 =4(r c sin 2 (20)+r A cos 2 (20)), 

72 = 4r c sin 2 (0)+r A . 

- classical- 



The measured probability consists of damped oscil- 
lations around a mean that approaches an asymptotic 
value. The decay of the oscillations is related only to the 
decoherence of p{t), while the relaxation of the popula- 
tions, encoded in Y(t), determines the slow time varia- 
tion of the mean. The dependence of the decay rates, 
7l; 72, on e is quite different from what we would have 
obtained by simply accounting for fluctuation of </(e), in 
which case 



7i = , 72 = 4r c sin 4 6 + T x sin 2 (26>) 
- classical energy fluctuations only - 



(8) 



This means that the effects of fluctuations of the state 
\S g ) cannot be neglected with respect to (exchange) en- 
ergy fluctuations. We also note that the asymptotic value 
of the probability is P(t) = 3/8, corresponding to a 
steady state density matrix with equally populated sin- 
glet states (\S g ), \S e )). This is a signature of the high 
temperature (fc^T 3> frui s ) limit, and is related to the as- 
sumption of classical gate voltage fluctuations. By con- 
trast, at low temperature we expect that only the lower 
singlet level is populated. In the experiment oscillations 
with J(e) ranging from tenths to few peV have been ob- 
served at ksT ~ lOpeV— . We therefore expect that the 
condition /cgT ~ Hlj s 3> J(e) is that obtained experi- 
mentally, in which case quantum state fluctuations are 
important and we explore their dependence on the tem- 
perature. 



IV. QUANTUM NOISE 



In order to extend our analysis beyond the high tem- 
perature limit we need to consider the quantum nature 
of the gate voltage fluctuations. We therefore modify our 
classical model (Eq. 



Here the operators V^fA) acting on the system are the 
same as in Eq. ([3]) while the classical fluctuators are re- 
placed by the operators A c{x) = £ 4 a c (\) t {b c (\)\ + b c(\) t ) 

of bosonic baths, -ffbathc(A) = Y,% ^i&c(A)]\(A)j- The ef- 
fect of the reservoirs on the dynamics of the electrons in 
the double dot is entirely characterized by their symmet- 
ric and antisymmetric spectral functions 



S c ± ) H = l/(27r) 



dte i " t ([A c(x) (t) > A cW (0)]±). (10) 



We take both baths to be at equilibrium at the same tem- 
perature ksT = 1/13. Their bosonic nature guarantees 
that S + (ui) = coth.(Ptkx)/2)S~ (u>), where S~(uj) is tem- 



(7) perature independent, S, x Juj) = J2i a c(A) 2 [<5(^ + Ui) 



S(uj - oji)}. 

We assume that the bath is weakly coupled to the sys- 
tem, and we determine the evolution of the density ma- 
trix to second order in a c (xy/u) s . Following the Bloch- 
Rcdficld approximation^ 1 ", we introduce a (short) bath 
correlation time, f, characterizing the typical time scale 
at which any correlation of the system and the reservoir 
disappears. The time evolution of the reduced density 
matrix of the system, p, coarse grained at time scales 
At ^> f, is Markovian. It is determined by the first order 
linear differential equation 



dtPa,, 



-ILUa.b Pa,b 



b,c,d Pc,de- 



(11) 



written in the basis of eigenstates of Ho where hui a ,b — 
(a\Ho\a) — (b\H \b) and TZ a ,b,c,d is the Bloch-Rcdficld ten- 
sor: 



Tia.,b,c.d = E y^^b4(a\Vj\ri)(n\Vj\c)gj(u; Ctn ) 

j=c,\ L n 

-(a\V ] \c)(d\V ] \b)g ] {^ a ) - {a\V 3 \c){d\V 3 \b) g] M 



■^Sa,c{d\Vj\n){n\Vj\b)gj{w n , a 



(12) 



H = H + V C A C + V\A\ + H] 



bath c 



bath A • 



(9) 



with 5 H = l/A[Sf(u))+Sj{u))]-i J WL Pdx/(2TT) (S+(x)+ 
Sj{x))/(x - oj). At times t > h/J(e) > H/w s , ne- 
glecting terms of order 0(J/u) s ), the sum in Eq. (TIT]) 
involves only terms such that h(u) a> b — uj c ,d) -C 0. Ex- 
plicitly the only relevant entries of the Bloch-Rcdfield 
tensor are: Ks g ,s g ,s g ,s g = -fcs e ,Se,S g ,S g , Rs e ,Se,S t ,s e = 

-^S s ,S 9 ,Se,S e , K Sg ,T L ,S g ,T L = T^T L ,S g ,T L S g - lt foll ° WS 

that state population and coherency evolve indepen- 
dently of each other and therefore the density matrix 
has the same form presented in Eq. ([5]), with different 
functions X(t) and Y(t). 

Once the expression of p(t) is known, it can be used to 
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calculate the survival probability. 

P qm (r) = I [3 + + tauh ^) (1 - e"^) 

+4 cos(( J(e)/ft + A, 7 )r)e-T 2r ] , (13) 

71 = 27r[sin 2 (260S'+(w s ) + cos 2 (26»)S'+(w s )]/ft 2 , (14) 

72 = 7T [sin 2 (20) (5+ (w.) - + cos 2 (20) 

- S£(u.)) + (cos(2#) - l) 2 5+(0) 

+ sin 2 (20)5+(O)]/(2n 2 ), (15) 

and n 2 Aj = [sin 2 (2<9) J R P oIuj/lu (S x (uj) - S~(u)) + 
J R PdLo/(u + uj s )(sm 2 (29)(S+(co) + 5~H) + 
cos 2 (26>)(S^H + 5^(w)))]/2. Aj is a shift in the 
frequency of the Rabi oscillations that can be neglected 
compared with J(e)/H, consistent with our second order 
perturbation expansion. The Bloch-Rcdficld approxima- 
tion employed implies that Eq. (|13[) is valid in the limit 
7i(2) ^ ^s,^/t. We note that, unlike 72, 71 depends 
only on the symmetric (classical) correlators, S^^(uj). 
72 consists of contributions from the bath correlation 
function at frequency oj s (which describes the relaxation 
process between the two singlet cigenstates and the 
corresponding contribution to the dephasing), and from 
the zero frequency correlation function (corresponding 
to the contribution of pure dephasing) 1 ^. 

In the high temperature limit, /3tkJ s — > 0, S~,y.(ui s ) 

is negligibly small as compared with S^^(uj s ), we thus 
expect a classical result. If we furthermore assume a 
Ohmic bath, i.e. S~,y.(u) = a c (\}h 2 u>, such that at high 

temperature S^Jw) ~ ^jm(O) for < w < lo s , Eq. (jT5|) 
reduces to Eq. ([6]) with 

71 = 2*/h(S+(0) sin 2 (20) + S+(0) cos 2 (20)) , 
72=V(2^)(45+(O)sin 2 (0)5+(O)). 

- quanta! high T - 

We in fact recover the classical result for white 
noise fluctuations (Eq. (O) by identifying T C ( A ) = 
7r / '(2H 2 )S't x JuJ s ). At low temperature, (3hw s ^> 1, 
S5y.(u> s ) ~ S~^ X j(ii! s ) and the quantum nature of the 
bath becomes important. The rate 71 disappears from 
the expression for P(t) and we also note the finite fre- 
quency contribution to 72 vanishes. Only the zero fre- 
quency component of the spectral density of the bath 
(which is responsible for pure dephasing) survives, 

72 = 7r/(2ft 2 )(4S'+(0) sin 4 + S+(0) sin 2 (2(9)) . (17) 

- quantal low T - 

In this limit the dependence of 72 on 9 can be explained 
entirely in terms of classical fluctuations of the oscillation 
frequency, J(e) (cf. Eq. ([5])). Indeed, the effects of fluc- 
tuations of the state \S g ) do involve transitions between 
the latter state and the singlet excited state, yet these 
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Figure 2: (color online) Surviving probability versus time r 
(Eq. (O) at J(e) = 3.9 ^eV (black solid line), J(e) = 2.4 /xeV 
(orange dashed line) and J(e) = 2 pieV (cyan dotted line). 
For all the curves J(0) = 5 peV, T = 135 meV as estimated 
from the experiment, and a c = a\ = 7 x 1CF 3 obtained by 
fitting with the experimental curves. The high and low tem- 
perature asymptotic values are depicted (dotted gray lines). 
Insets: Q = Jji^l-Khrj-z) oc (number of visible oscillations) as 
a function of J(e) for three possible scenarios: (a) a c = 0.02, 
a x = 0; (b) a c = q a = 7 x 10 -3 ; (c) a c = 0, a x = 7 x 10~ 3 . 
Different curves correspond to different temperatures accord- 
ing to their labels in the inset (c); in all plots J(0) = 5 /j,eV. 

transitions are exponentially suppressed by e~ l3uJs . Note 
however that experimentally the regime (3ftw s < 1 can be 
reached, hence the effect of fluctuations of the state \S g ) 
can be of interest. In particular this fluctuations affect 
the steady state value of the survival probability which is 
a function of (3u> s , Pqm( T ~ ¥ 00 ) = l/8[3 + tanh(/37ia; s /2)] 
(cf. Fig. [2]) and can be directly observed in the experi- 
ments. 

A comparison with the experimental results of Ref. [H 
is obtained assuming Ohmic baths with spectral den- 
sities S~,y,(u>) — a c (\)h 2 uj. These properly describe 
charge fluctuations due to the external circuit controlling 
the gate voltages. An analysis of the possible scenarios 
etc p oe\ , shows that the experimental fact that the num- 
ber of visible oscillations as function of J(e) is constant 
is correctly reproduced for a c < a\ (cf. insets in Fig. 
The time dependence of P qm (T) is depicted in Fig. [2] for 
different values of J(e). We obtain a fit with the experi- 
mental data for a c = at\ = 7 x 10~ 3 , consistent with the 
strength of electromagnetic environment in other solid 
state systems^. 

V. EXTENDED MODEL 

The previous analysis is now extended to include the 
lowest energy triplet state in the charge configuration 
(0,2), \Tq). This might be necessary if the energy of 
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|Tq) is comparable with the Coulomb energy Eq. The 
new Hamiltonian reads, H' = H + E (X t \T ) (Iq| + 
h.c.) + Eq(S — X s e) |Tq) (Tq|, where S is the excitation 
energy to the triplet state in the (0, 2) configuration. 
Owing to electron tunneling, At, the two triplet states, 
\T Q ) and \T£), do hybridize (cf. Fig. [Hc))±i. The en- 
ergy spectrum of singlet and triplet states in the subspace 
S z = is depicted in Fig. 0Jc). The energies of the hy- 
bridized triplet states, \T g (e)) = — sui^ITq) + cos ip \T ) 
and \T e (ej) = cos^ |T,5) + sin^ \T ), are E 5 + = 
Eq5 + Eoy / (e — S) 2 + A 2 ] respectively, with tan 99 = 
(\/(e - 5) 2 + Xf + e - S)/X t . The Hamiltonian for this 
model includes four parameters, e, 5, A s , At which depend 
on three fluctuating gate voltages only, Vl,Vr,Vt- In 
principle it is possible to determine a correlation ma- 
trix for the fluctuations of the parameters in the Hamil- 
tonian by considering a specific potential form for the 
double dot. Instead we assume that Vr affects only 
the tunneling matrix elements, which is reasonable for 
weak tunneling. The fluctuations of At — * At + £t(t) 
and A s — > A s + £,s{t) will then depend on the same 
bath and will therefore be correlated, £ t (i) = f£ s (t) 
with / = (dX t /dV T )(dV T /dX s ). At the same time the 
gate voltage difference Vl — Vr will affects only e. The 
density matrix now evolves as (cf. Eq. ©) p'(t) = 
p(t) + l/2W(t) \T g ) {T g \ + 1/2(1- W(t)) \T e ) (T e |, and the 
probability of finding the system in the |||) at time r is 
P' qm {r) = P qm {r) - 1/8(1 - tanh(/?^ t /2))(l - e"^), 
with 73 = 2ir/h 2 [sm 2 {2ip)S+(uj t ) + P cos 2 {2<p)S+ (uj t )} 
and 72 replaced by 7 2 = 7i/4(l — tanh(/3fiw s )) + 73/4(1 — 
tanh(/3fiw t )) + 7r/(2ft 2 )[(cos 2 (26>) - cos 2 (2<p))S + (0) + 
(sin 2 (26>) - /sin 2 (2^))S^(0)]. The physical mechanism 
that induces decoherence in P(r) is the same described 
in the previous paragraphs. Similarly to the decoher- 
ence in the singlet subspace, fluctuations in \T g ) do in- 
volve now the exited triplet state |T e ), an effect that is 



small in e ~P h ^tl 2 . Note that, even at low temperature, 
fcfiT <C tkjJt, fluctuations in the energy J(e) (cf. Fig.QTc)) 
modify P qm (r). Remarkably in this case, even in the 
presence of a "sweet point" (9 = <p)^, while d e J = 0, 
fluctuations of the tunneling rates are important due to 
the difference between electron tunneling in the triplet 
and singlet states. 



VI. CONCLUSIONS 

We have presented here a simple model describing the 
effect of charge fluctuations on Rabi oscillations between 
spin states |t!) and |||) °f electrons in a double dot. We 
have accounted for decoherence effects due to both energy 
and quantum state fluctuations, by including the quan- 
tum effects of a fluctuating environment within the Born- 
Markov approximation — Eqs. (fT51[T5|) . We have shown 
that not only in the high temperature limit does the 
result reproduce that of classical fluctuations (compare 
Eq. (fl6|) to Eq. (J7J)), but also the low temperature result 
has a classical interpretation in terms of energy fluctu- 
ations only (not state fluctuations, compare Eq. (jT5J) to 
Eq. ©). In fact the role of the state fluctuations is sig- 
nificant at a temperature that exceeds the singlet excita- 
tion energy, a regime which is accessible experimentally. 
Note that at high temperature the "classical limit" refers 
to classical environment induced fluctuations. The latter 
can still cause fluctuations in the quantum state of the 
system. At low temperature state fluctuations are frozen 
out. 
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Israel BSF, the ISF of the Israel Academy of Sciences and 
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